Graphene plasmons and retardation: strong light-matter coupling 
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Abstract -We study the retardation regime of doped graphene plasmons, given by the nominal 
crossing of the unretarded plasmon and light-cone. In addition to modifications in the plasmon 
dispersion relation, retardation implies strong coupling between propagating light and matter, 
even for homogeneous graphene, which opens up the possibility of efficient plasmonics in simple 
graphene devices. We exemplify this enhancement in a double-layer configuration that exhibits 
perfect (if lossless) light transmissions across a classically forbidden region, providing a simpler 
analog of the corresponding phenomenon in perforated metal sheets. We also show that (broad) 
Fabry-Perot resonances present without graphene turn into sharply peaked, quasi-discrete modes 
in the presence of graphene where graphene's response function is given by the typical Fano 
lineshape. 



Introduction. — In addition to its amazing transport 
and mechanical properties pQ, graphene's optical behavior 
is also notable. Absorption, for instance, has the exper- 
imentally observed [2j [3], universal value « ira for light 
in the visible spectrum, depending on the fine structure 
constant a, but not on material's properties. Several pro- 
posals and/or realizations highlight graphene potential in 
optical and communication technologies, including resis- 
tive touchscreens of transparent and flexible displays, [4]l5] , 
opto-electronic devices such as photodetector [6], broad- 
band absorber |7], mode- locked laser due to current sat- 
uration [8], or as a highly efficient fluorescence quencher 

[MS]. 

Doped graphene has attracted much attention recently 
as a suitable candidate for noble metal's replacement in 
the active field of plasmonics [T6HT9] . In addition to their 
intrinsic importance as a probe for dynamics, graphene 
plasmons [20H22] (GPs) offer a number of advantages such 
as frequency tunability, long life-times, and spatial con- 
finement due shorter (than light) wavelengths. They have 
been observed via electron energy loss spectroscopy [23]l24] 
and near-field nanoscopy. [25H27] 

Unfortunately, GPs do not couple easily to propagat- 
ing electromagnetic modes over most of its range due to 
the large momentum mismatch, complicating their ma- 
nipulation. This problem could be overcome by breaking 
the conservation of parallel momentum either with con- 



fined geometries [19JGE] or artificial periodicities 
However, there is a regime where even homogeneous GPs 
must couple strongly to (propagating) light: the retarda- 
tion limit. It applies for frequencies below a characteristic 
crossover scale, which can be taken as the crossing between 
the nominal, unretarded GP dispersion and the light-cone. 
There, phenomena associated with strong light-graphene 
coupling take place. 

Here, we demonstrate this enhanced light-graphene cou- 
pling and apply it to a double-layer graphene |32 ? 33] ar- 
rangement possessing extraordinary transmission. This 
term was originally coined to describe the enormous trans- 
mission experimentally observed through periodically per- 
forated metal sheets [34H36] , a situation where the naive 
expectation would have prescribed just the opposite. An 
explanation was provided in terms of the excitation of sur- 
face plasmons, with the result of enhanced transmission 
(perfect in the absence of absorption) through a nominally 
opaque region. In our case, the resonant coherent excita- 
tions of the graphene layers allow the enhanced transmis- 
sion of photons through the central, classically forbidden 
region for photons, in direct analogy with the metallic 
case. In the propagating region, the strongly quenched 
transmission turns into perfect transmission at sharply 
peaked resonances related to the original Fabry-Perot res- 
onances, where graphene's response function is given by 
the typical Fano lineshape due to the nearly discrete na- 
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ture of the spectrum in the central slab. 

The Letter is organized as follows. We first study 
the retardation limit and its effect on the GP dispersion 
for single-layer graphene. Then we consider the retarda- 
tion regime for double-layer graphene separating different 
dielectric media where extraordinary transmission takes 
place. Finally we summarize our results and an appendix 
details the GPs dispersion relation for the double-layer 
arrangement in the retardation limit. 

Retardation limit. Single-layer graphene. — The 

standard expression for a 2d plasmon, uj p oc y/q, assumes 
instantaneous Coulomb coupling between charges. [201 121] 
Therefore, it cannot be correct when the nominal plas- 
mon dispersion meets the light-cone, uj p < cq. This is the 
regime where retardation effects matter and strong light- 
graphene interaction takes place, as we now show. 

Consider a graphene layer sandwiched between two di- 
electrics with permittivities ei^ = e lj2 £o anc ^ e i > e 2- 
The longitudinal current response to an (in-plane) longi- 
tudinal, external vector potential is given by the standard 
RPA expression 
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Fig. 1: Left panel: Longitudinal plasmon exact dispersion re- 
lation (black continuous line), compared to the instantaneous 
approximation (blue dashed line) for doped graphene between 
medium 1 (ei = 5) and medium 2 (e2 = 1), light-cones also 
shown. Right panel: As in left panel for the double-layer 
graphene showing the in-phase (upper curve, black) and out- 
of-phase (lower curve, magenta) plasmons for 63 = e± and layer 
separation z = z c /2. Inset: schematic geometry and transmis- 
sion setup. 



where xi — Gjj is the retarded Green function for the for doping level n ~ 10 12 cm 
longitudinal current, x° its non-interacting (bare) ver- 
sion, and d\ = Qaa that of the (in-plane) longitudi- 
nal vector potential (in the absence of graphene). For 
graphene, we use the well-known results [2Q | [2T | [37] for x°, 
whereas for the photon field, one straightforwardly finds 
QW2 c 
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GPs are the poles of Eq. [TJ leading to the following 



(implicit) dispersion relation, including retardation 
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In the unretarded limit, c — >• 00, this gives the known 
square-root dispersion. In contrast, the exact solution re- 
places this behavior with a linear dispersion which merges 
with the slower medium (1) light-cone below a characteris- 
tic crossover frequency, as illustrated in Fig. [I] (left panel), 
and summarized as follows: 
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where uop^F and v are graphene's Fermi frequency, Fermi 
momentum and Fermi velocity, c\ = c/y^el is medium 1 
(slower) light velocity, and a g = represents graphene 
fine structure constant. The crossover between regimes 
takes place for frequencies which roughly corresponds to 
the intersection of the unretarded GP and light-cone dis- 
persions. This scale is given by uo c ~ auop for reasonable 
dielectric constants. Typical frequencies are v c ~ 200 GHz 



v c ~ 600 GHz for doping 
level n ~ 10 13 cm -2 , reaching the technologically impor- 
tant THz regime for n ~ 10 14 cm -2 . 

Therefore, for uj < uj c retardation always matters. Not 
surprisingly, this linear regime is also the region of strong 
graphene-light coupling. This is immediately seen by look- 
ing at the reflection and transmission amplitudes for the 
(in-plane) longitudinal vector potential upon passing from 
medium i to j, given by 



l+Ti 



(4) 



For interband transitions at frequencies above 2cjf, 
graphene terms in Eq. [4] are minute, leading to the univer- 
sal 2.3% weak absorption in vacuum, for instance. On the 
other hand, for frequencies below uo c1 graphene response 
starts to dominate in Eq. [4] implying strong radiation- 
graphene coupling. For instance, in the limit u <C uo c the 
reflection amplitude becomes 



r = -l + 0( — ). 



(5) 



r = — 1 implies perfect reflection or, equivalent ly, a zero 
photon field at graphene position as a boundary condi- 
tion. Therefore, the (small) parameter — ~ — — of our 
strong coupling regime also measures the departure from 
this zero-field boundary condition. 

Double-layer graphene. Extraordinary transmis- 
sion. — We consider two identically doped graphene 
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sheets, separating three dielectrics with permittivities 
e i,2,3, see Fig. [T] (right panel), and impose left-right sym- 
metry choosing e 2 < e\ =63. Under these conditions, for 
large enough incident angles the central region does not 
support propagating modes: light must tunnel through 
this evanescent region. Doped graphene enhanced re- 
sponse in the retardation regime can make this tunneling 
perfect, as we will show. 

Reflection and transmission coefficients can be obtained 
from the corresponding amplitudes, written as 



n,3 = r 12 + 



£i2r 2 3^2ie 2q * z 

-2q' 2 z 



1 - r 2 ir 23 e 
ti2t2se~ q ' 2Z 



(6) 



c l,3 



1 - r 2 ir 23 e 



-2q' 2 z 



for light incoming from medium 1 and being transmitted 
to medium 3 (see Fig. [T]inset), with tij and rij taken from 
Eq. H 

In order to display the connection between light 
propagation and graphene behavior, we also calculate 
the double-layer (longitudinal) current-current response, 
given by the following matrix generalization of Eq. [T] 



Xl = 0--e 2 x°di)- 1 x°, 



(7) 



with the 2x2 matrix x° — diag(xi,X 2 ) representing the 
non-interacting (longitudinal) response of graphene layers 
1 and 2, whereas the photonic matrix d = (dij) is given 
by dn = d\ (1 + ff >3 ), d 2 2 = d 3 (1 + fg^), and d i2 = 

di?f >3 = ^3 ^3,1 = d 2i- Here, d x = d 3 = and 

f°j (t°j ) correspond to the expressions of Eq. [6] evaluated 
for the dielectric geometry of Fig. [TJ but without graphene 
layers (x° = in Eq. 2]). For our left-right symmetric 
arrangement, the in-phase and out-of-phase components 
of graphene currents, j± = ji ± j 2 , diagonalize with 
corresponding (complex) eigenvalues that we denote by 
X++ and X---) respectively. 

Genuine (non radiative) GPs correspond to the poles 
of xi m the totally evanescent regime, that is, the region 
to the right of the slower light-cone in Fig. [TJ In the 
Appendix we show that there always are two GPs, one 
for each component of the response. The in-phase GP 
(X++) merges with the corresponding limiting light-cone, 
whereas the out-of-phase GP (X--) either merges with the 
light-cone or develops its own linear dispersion at slower 
velocity if the separation between layers is below a critical 
value z c = z c kp, with z c = 62 tt— . Both are shown 

L L L ' L ei — €2 2 V OL 

in Fig. [D (right panel) where we have used z = z c /2, 
and permittivities e 2 = 1 (vacuum) and e\ = 63 = 5, a 
representative value for common substrates. For bulk BN, 
we have, e.g., 61,3 ~ 5, whereas 61,3 ~ 4 for few layer BN 
or Si0 2 . 

Though GPs remain in the evanescent region of both 
media, the presence of graphene also strongly affects prop- 
agating light, leading to resonances with perfect transmis- 
sion. At first glance, this statement seems incompatible 
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Fig. 2: Top panel: Transmission as a function of incoming 
angle for frequency uj — (a/2)ujF and double-layer graphene 
separation zkF = 20c/(va) for permittivities 61,3 = 5, €2 = 1. 
Dashed (black) line: same as before in the absence of graphene. 
Bottom panel: Spectral density of the in-phase (solid, blue 
line) and out-of-phase (dashed, red line) double- layer graphene 
response. 



with the (almost) perfect reflective character of a graphene 
layer in our strong coupling regime. Yet, the zero field 
boundary condition also decouples the photon field in the 
central slab from the propagating modes outside, leading 
to the standard discrete spectrum for a particle in a box. 
Now, the small parameter ^- ~ measuring the depar- 
ture from perfect reflection is also a measure of the (small) 
coupling between the discrete modes of the central slab 
and the propagating modes outside, leading to (radiative) 
resonances. Drastic changes are expected close to these 
resonances and, indeed, we will see that the otherwise ex- 
pected strong reflection turns into perfect transmission. 

These resonances are shown in Fig. [2] (top panel), 
where we present the transmission for layer's separation 
zkp = 20c/(va) and frequency uj = (ql/2)ujf- For a 
doping level of n ~ 10 12 cm -2 , this would imply a fre- 
quency v ~ 100 GHz, and separation z~5x 10 _3 m. For 
n ~ 10 14 cm -2 , we have v ~ 1 THz, and z~5x 10 _4 m. 
Although the existence of perfect transmission at Fabry- 
Perot-like resonances is generic for this geometry even 
without graphene (see dotted line in top panel of Fig|2]), 
these resonances are sharply modified by its presence. The 
otherwise strongly reflective graphene confines the trans- 
mission to narrow resonances close to the region where 
the almost discrete central slab modes are excited. The 
narrow width of these resonances is a measure of the long 
lifetime of the central photonic modes, weakly decaying 
into the radiative modes outside. A direct evidence of the 
role of these (almost) discrete spectrum is provided by the 
spectral densities of the in-phase (x++) an d out-of-phase 
(X--) graphene response, also shown in Fig. [2] (bottom 
panel). Up to an overall scale, these spectral densities co- 
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incide with the spectral densities of the photonic field and, 
therefore, the sharp resonances in graphene match the ex- 
pected sharp spectral densities of weakly decaying, almost 
discrete photonic modes. Notice the one-to-one correspon- 
dence between response and transmission peaks, empha- 
sizing the role of these modes in overturning the otherwise 
generic strong reflection of graphene in the strong coupling 
limit. 

The spectral functions of the bottom panel of Fig. [2] 
exhibit a marked asymmetric Fano-like profile, vanishing 
precisely where grapheneless transmission is one. This 
Fano-like profile can be expected, given the key role played 
by the almost discrete modes of the central slab. It can 
be quantitatively explained noticing that 



1 



X±± 



l-e 2 x^ii±d 12 )' 



(8) 



but the photon propagator in the absence of graphene van- 
ishes at the corresponding grapheneless resonances as 



du ± d±2 ~ a(Ss) + ib(Ss) 2 



(9) 



where a and b are real constants, and Ss is the depar- 
ture of sin(0i) from the value corresponding to the dis- 
crete modes upon strictly enforcing zero-field boundary 
conditions, given by exp(— 2q' 2 z) = 1. Notice that this 
zero-field boundary condition for discrete modes in the 
presence of graphene coincides with the perfect transmis- 
sion condition in the absence of graphene, explaining the 
vicinity of both points. The small shift of the actual solu- 
tion being the consequence of the small leaking of dis- 
crete modes into the outside continuum. This leaking 
also accounts for the small width of the resonance, and 
leads to a characteristic resonance asymmetry, given by 

~ Im X±± ~ (Ss-sly+i^ss)^ where s * ~ is the res " 

onance shift from the nominal discrete modes graphene, 
playing the role of small parameter in the strong coupling 
regime, and I being (up to a factor of order one) the central 
slab width in units of the minimum required for the emer- 
gence of the first nominal discrete mode. This expression 
can be approximately recast in the more common Fano 
lineshape [38] , 



/\2 
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(5s>y + ( 7 /2)2 



(10) 



with the identifications Ss' = 5s — s*, the width 7 ~ /s 2 , 
and the Fano asymmetry parameter Q = ^ with value 
Q ~ 2 — 3, for fig. [2] (bottom panel). 

Interestingly, the rightmost resonance in Fig. [2] oc- 
curs for sin 6^ > ~ 0.447, that is, in the evanes- 
cent regime for medium 2. In fact, we can tune the 
arrangement to make this last resonance the only one 
present. This is shown in Fig. [3j where we present 
the transmission calculated with the same frequency as 
before uj = (a/2)cjp, but for three shorter separations 
given by z\kp = 2c/(va), z<i = zi/5, and Z3 = zi/10 




sin 6. 



Fig. 3: Transmission as a function of incoming angle for fre- 
quency uj = (a/2)u)F, permittivities 61,3 = 5, €2 = 1, and 
double-layer graphene separations: ziJzf — 2c/(va) (black), 
Z2 = zi/5 (red), and Z3 = zi/10 (blue). The vertical dashed 
line marks the onset of the evanescent regime in the central 
slab. 



(v ~ 100 GHz and z x ~ 5 x 10- 4 m for n ~ 10 12 cm- 2 ). 
One can observe a perfect transmission peak within the 
classically forbidden region for photons. This resonance is 
rather sharp close to the propagating boundary, moving 
deeper into the evanescent region and broadening with de- 
creasing double-layer separation. Notice that, for the cases 
shown in Fig. [3j the transparency window is virtually con- 
fined to the evanescent regime. Although a finite photon 
tunneling always exists, the fact that the transmission be- 
comes perfect in the evanescent region [39] is entirely due 
to the presence of doped graphene sheets, and is associ- 
ated with a resonance in their out-of-phase response (see 
bottom panel of Fig. [2]). This perfect transmission peak 
within the evanescent region extends all the way down 
to zero frequency, as shown in the top left panel of Fig. 
[4] For higher frequencies, this resonance approaches the 
upper light-cone and merges with the lowest propagat- 
ing Fabry-Perot like resonance, as seen in the bottom left 
panel of Fig. [4j 

Our aim so far has been to expose the strong light- 
graphene coupling of the retardation regime in a sim- 
ple manner. Therefore, we have ignored absorption. 
Graphene losses will degrade the perfect nature of trans- 
mission, making the possible observation of coherent ef- 
fects more likely in reflection rather than transmission. 
Although we are not aware of any intrinsic limitation in 
the minimization of losses attainable in doped graphene, 
the low frequencies of the strong-coupling regime might 
pose a severe handicap. Nevertheless, absorption opens 
new possibilities: the sharp response of graphene in the ab- 
sence of absorption suggests a correspondingly enhanced 
absorption when losses are allowed. This is the case as 
the upper right panel of Fig. [4] shows, where the ab- 



p-4 



Graphene plasmons and retardation 




Fig. 4: Top left panel: Transmission in the uo — q plane showing 
the perfect resonance within the evanescent region for double- 
layer separation corresponding to zik,F — 2c/(5v a). Bottom 
left panel: As in top left panel for an extended region of the 
uj — q plane to include additional Fabry- Perot- like resonances. 
Top right panel: Absorption when graphene losses are included 
with ujt = 0.65 for the same double-layer separations (and 
colors) as in Fig. [3j Bottom right panel: Absorption for fixed 
double-layer separation z\k,F — 2c/(va) and graphene losses 
lot = 2 (black), ujt = 0.5 (red), and ujt = 0.05 (blue). 

sorption for a lifetime r ~ 10 _12 s, present in suspended 
graphene and corresponding to ujt = 0.65, is plotted for 
the same double-layer separations of Fig. [3j Indeed, one 
observes an overall high absorption, reaching values above 
90% precisely in the evanescent region. From this perspec- 
tive, our arrangement is closely connected to that of Ref. 
[40] , providing a complementary and potentially simpler 
alternative to absorption enhancement based on periodic 
patterning [29H5T] . It is similar, though, to the previously 
suggested enhanced absorption of graphene placed in a 
(double) Fabry-Perot cavity. p[TlH2] 

In the bottom right panel of Fig. 2J the absorp- 
tion for different intrinsic damping rates corresponding 
to lot = 2,0.5,0.05 is shown where the last value would 
assume a phenomenological relaxation rate 7 of the po- 
larization with h'y ~ 10meV, present in graphene on a 
Si02-substrate. Notice that the double-layer absorption 
can display a non-monotonous behavior with the intrinsic 
single-layer graphene losses, providing a richer degree of 
control. 

Summary. — We have studied the retardation regime 
of doped graphene plasmons, given by (uo/uop) ^ ol. Apart 
from modifying the unretarded GP dispersion behavior 
from ^fq to such a limit marks the onset of strong 
radiation-graphene coupling and we have exhibited such 
enhancement in the simplest double-layer arrangement, 
i) We first observed that the transmission is strongly 
quenched except at sharply peaked resonances, consistent 
with the nearly discrete spectrum of the strongly confined 
field in the central slab brought about by graphene, whose 



response function is given by the typical Fano lineshape. 
ii) We also discussed that the transmission, perfect if losses 
are ignored, can be confined to the classically forbidden 
region for photons between graphene layers, providing a 
direct analog of the extraordinary transmission through 
perforated metals in a conceptually simpler setup, iii) 
Finally, we showed that enhanced transmission becomes 
increased absorption when losses are allowed, displaying 
non-monotonous behavior. Graphene's doping tunability 
in this or similar setups should thus open up new ways to 
efficiently control the flow and/or absorption of radiation. 

* * * 

This work has been supported by FCT via grant 
PTDC/FIS/101434/2008 and by MICINN via grant 
FIS2010-21883-C02-02. 

Appendix. — We outline here the calculation of 
double-layer GPs in the retarded regime. We take 63 = 
ei > 62- Therefore, genuine GPs are the response poles 
located to the right of the slower (outer media) light-cone. 
These are given by the solutions of 

l T r 21 e-^ = 0, (11) 

corresponding to the in-phase (x++) an d out-of-phase 
(X— ) response components, respectively. Taking ri\ from 
Eq. [4] deep into the retarded regime, one finds the in-phase 
GP as 

ei^^ 4 7T^xVi, (12) 
where we have used dimensionless magnitudes: Co = 
uo/uof-, otg = ca, c = c/v, cf = c 2 /ei, q = q/kF, 
q[ = \fq 2 — 6j 2 /c 2 , with single-layer graphene response, 
X° = 7T _1 , taken in the local approximation, excellent for 
the retarded regime. In this limit, oj + coincides with the 
GP pole of an isolated graphene layer in the boundary be- 
tween two (semi-infinite) media 1 and 2. Therefore, one 
easily sees that it merges asymptotically with the outer 
media light-cone, as in Eq. [3j 

0J + = c 1 q + O(f). (13) 

Analyzing the plus Eq. [11] for the out-of-phase GP, a 
critical layer separation emerges, z c , given by 



that separates the following two regimes 

Co. = ( C >1+ ! > * c , (15) 

\ v_q + 0(q 3 ), z< z c v ' 

where 

v 2 = C ^—~ . (16) 

" i + (i-a)(^-i) 

We conclude that the out-of-phase GP develops a linear 
dispersion relation that either merges with the outer media 
light-cone (and, therefore, with Co + from below) for z > z c , 
or shows a lower velocity v_ < c\ for z < z c , as stated in 
the paper. 
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